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ABSTRACT 

A statistical model of spatial context 
is described and procedures for classifying 
remote sensing data using a context clMsi- 
fier are outlined. Experimental results 
aire presented. Because the con^utational 
requirements of the context classifier are 
very large, its implementation on peurallel/ 
pipelined multiprocessor systems is being 
investigated. Some of the special consider- 
ations necessary for such implementations 
are described, with particular reference to 
implementation on an array of Control Data 
Corporation Flexible Processors. 


I. INTRODUCTION 

For more than a decade, efforts to 
extract information from multispectral 
remote sensing image data have proved 
increasingly successful. To a large extent, 
these efforts have focused on the applica- 
tion of pattern recognition techniques to 
the multispectral measurements made on 
individual ground resolution elements; l.e., 
scenes have been classified plxel-by-plxel 
based on the measurement vectors associated 
with the individual pixels ' . Progress has 
been achieved through development of 
increasingly sophisticated methods for ex- 
tracting information from the spectral do- 
main to characterize the classes of inter- 
est. 

However, there are many applications 
for which the classes of Interest can be 
better characterized if the spatial infor- 
mation in the remote sensing data is 
utilized in addition to the spectral 
measurements . Characteristic spatial 
features include, for example, shape, 
texture, and structural relationships. 


This work was sponsored in part by the 
National Aeronautics and Space Administra- 
tion under Contract NAS9-15466. 


Some interesting and useful research has 
been accomplished in recent years in the 
direction of Incorporating spatial informa- 
tion into the data analysis process^'*'*. 

One way to approach spatial information 
in image data is to recognize that the 
ground cover associated with a given pixel, 
i.e., its "class,” is not independent of 
the classes of its neighboring pixels. 

Stated in terms of a statistical classifica 
tion framework, we may have a better chance 
of correctly classifying a given pixel if 
we take account of not only the spectral 
measurements associated with the pixel 
itself but of the measurements and/or 
classifications of its "neighbors" as well. 
Notice that at some point we must make 
clear how "neighborhood" is to be defined. 

If the objects in the scene tend to be 
rather large relative to the resolution of 
the sensor, i.e., each object is likely to 
consist of many spectrally similar pixels, 
this fact can be exploited nicely by apply- 
ing a combination of scene segmentation 
techniques and sample classification (somer- 
times called "per-field" classification)’.: 
More generally, the image can be considered 
a two-dimensional random process and the 
characteristics of this process incorpora- 
ted into the classification strategy. This 
is the objective of the approach described 
here, in which a form of compound decision 
theory is employed to improve scene classi- 
fication through use of a statistical 
characterization of context. Our work is 
an extension of an idea by Welch et. al.’ 
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As increasingly complex forms of data 
and data analysis methods are employed, 
the computational requirements tend to 
become more demanding. Although improve- 
ment in the raw speed of digital computer 
components can be exploited to some extent 
to meet these requirements, it is clear 
that evolving computer architectures, 
especially those Involving multiple proc- 
essing elements, have much to offer. The 
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cimtaxt cl«tfifi«r d«*erlb«d have hai com- 
putational raquiramanta which art aavara and 
bacoma mora ao aa tha aisa ot tha contaxtual 
naighborhood ia axpandad. It ia a natural 
candidatai tharafora* for multiprocaaaor 
implamantation. 


IX* THE CONmt CLAIRXFXBit 

Tha Imaga data to bo claaalflad in 
aaaunad to ba a two-dimanaional N| x Nj 
array of multivariata pixala. Aaaociatad 
with tha pixal at "row i" and "column j" ia 
tha multivariata maaauromant vactor X^j r R** 
and tha trua atato or olaaa of tha pixal 
e ^ ■ (wx r • • • I Um) • Tha maaauramanta hava 
clasa-conditional danaitiaa p(X|u|,)* i * 1* 
2, m* and aro aaaumad to ba claaa- 

condltionally indapandont. Tha objactiva 
ia to claaaify tha N ■ x N 2 obaarvationa 
in tha array. 

■nia action (olaaaification) datarminad 
by tha claaaifiar for pixal (i* j) ia 
danotad by a^ j e 0. To puraua a Bayaaian 
(minimum riaxlf atratagy , lat tha loaa 
incurrad by taking action aj,j whan tha trua 
claaa la ba danotad by a^j) for 

aoma flxad fton-nagativa function !.(•, •). 

Tha avarage loaa incurrad ovor tha N claaai- 
ficationa in the array ia defined to be 

I *ij^ • 

in the moat general case, the action a^j 
may depend on all of the ohaervationa in 
the array. Let X denote thia "vector of 
vectors" » then the expected loss is 


S MOij, .„UII 

( 2 ) 

• - I I, aij<5))3 

1 # j 

and we would like to have a decision rule 
(the rule of choosing a^j based on X) 
which minimizee R(X) . Note that the expec- 
tation is with respect to 6ij. 

When context is ignored, the action 
(classification) depends only on the 
measurement vector Xij of the pixel to be 
classified, in which case axj (X) ■ aij(Xxj). 
For our present purposes, however, we'^wanr 
to incorporate some neighborhood information 
in the decision process, so we define a 
neighborhood, the "context," consisting of 
an arrangement of p pixels such as shown in 
Figure 1. The arrangement actually used 
will be based on physical and other practi- 
cal considerations related to the environ- 


ment and application. Lat Xij ba a p-vaotor 
of meaaurament vaotora aasooiatad with 
pixal (i, j) to ba olaaaifiad and lat Sm 
ba tiia oorraaponding p-veotor of aotual 
olaaaaa. Tha function aij(jy,j) mapa 
p-vactora of obaarvationa into single 
olassaa (i.a,, olaasifiaa pixal (i, j) 
based on X^j). Tha axpaotad loss over tha 
full array is 

R(X) - I •ij<3Sij>>3 (3) 

Furthermore, if L(«, •) is taken to ba tha 
0-1 loss function (no loss for a correct 
classification, unit loss for an error) and 
tha maasuremants in a neighborhood ara 
assumed to be class-conditional ly independ- 
ent, it can ba shown that Eg. (3) will ba 
minimized if every pixel is classified 
(action a is selected) so as to maximiza 

I r n f(X,|0,) 1 0**(ii^) (4) 
ix j e fiP, L J ^ 

®ij - • 

where Oj, and X, ara tha class and measure- 
ment of the itn pixal in tha p-array (in 
any convaniant ordar) , f(X|0) is tha class- 
conditional density of X, and GP(0£j) - 
6F(6i, @2' "" ^p) ' which ideally must be 
known for the type of scene to be classi- 
fied, but in practice must usually be 
estimated from an accurately classified 
sample of the scene or from an analogous 
scene of known classification. 

An experiment was formulated to 
investigate the extent to which this 
classifier model can utilize contextual 
Information in satellite-gathered remote 
sensing data. In order to avoid confounding 
other effects with the impact of context, 
it was decided to use a simulated data set 
generated as follows. A classification of 
multispectral remote sensing data was 
selected which had been judged to be very 
accurate (typically, produced by careful 
analysis and refinement of multitemporal 
data) . Such a classification could be 
expected to embody the contextual content 
of an actual ground scene. Based on the 
classification map and using the associated 
statistics of the classes (developed in 
producing the classification) data vectors 
were produced by a Gaussian random number 
generator and composed into a new data set. 
Thru the new data set had the following 
characteristics s 

1. Each pixel in the simulated 
data set represented the 
same class as in the "template” 


R(X) - E 
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olMBifieation. Th« taaplat* 
oottld b« oonsldarad th« Aground 
truth* tor th« now data aat. 

2. MI. elaisaa in tha data aat 
wara knowi and rapraaantad. 

3. All elaiaaa had multivariata 
Oauaalan diatributiona with 
atatiatioa typioal of thoaa 
found in raal data. 

4. All pixala wara olaaa<>oondi» 
tionally indapandant of 
adjaoant pixala. 

5. Thara wara no nixtura pixala. 

Although tha aimulabad data aonundiat 
of an Idaalixation of "raal* romo^'j aanaing. 
ita apatial organiaation ia oonaia^^tint with 
a raal world aoana and ita ovarall eharaetar- 
iatioa ara oonaiatant with tha oontaxt nodal 
aat out abova. Zn aaaanoa. than, what tha 
axparimantal raaulta baaad on tha aimulatad 
datA ahow ia tha affaotivanaaa of tha oon« 
taxt olaaaifiar glvan that tha undarlying 
aaaumptiona ara raaaonabla. Purthar axpari* 
manta ara raguirad to ganaraliaa tha oonelu- 
aiona of thaaa raaulta to raal data. 

Thraa data aata wara aalaotad to rapra* 
aant a variaty of ground oovar typaa and 
taxturaa. Data aat 1 ia agrioultural 
(Hilliaton. North Dakota) « with ground 
raaolution and apaet.ral banda appro: :lma ting 
thoaa of tha projaotad Landaat D Thamatio 
Mappar. Data aat 2a ia Landaat 1 data from 
an urban araa (Grand Raplda. Michigan). 

Data aat 2b ia from tha aama Landaat frame 
aa 2a. but from a loeala having aignifi- 
cantly diffarant apatial organiaation. Each 
data aat la aguara. 50 pixala on a aida. 

Flgura 2 ahowa tha achiaved olaaaifi- 
cation raaulta. Tha "no oontaxt" clanalfi-* 
cation accuracy ia plotted coincident with 
the vertical axia of each graph. Data aet 
1 waa claaalfied ualng aucceBaivaly 2. 4. 6 
and 8 neighboring pixala; data aata 2a and 
2b were claaaified uaing 2. 4. and 8 
neighboring pixela. The reaulta apeak for 
themaelvaa. The accuracy improvement 
raaulting from the uaa of contextual infor- 
mation la quite algnif leant. 

For thia experiment, the context dia- 
tribution GPce^j) waa aimply tabulated from 
the "template'” clasalfieatlon. But in a 
real data situation, such a template la not 
available (else there would be no need to 
perform any further classification) . one 
can envision a number of ways in which the 
p-vector distribution might be estimated for 
a remote sensing application. For example, 
it could be extracted from a classification 
of the same araa obtained previously. This 


would raquira that tha araa not have 
ehangad too graatly in ita olaaa maka-iqp 
slnoa tha aarliar data wara oollaotad and 
that tha aarliar olaaelfieation was raaaon- 
ably aeeurata. Or. the distributimt night 
ba <^tainad fren a olaasifloation of any 
■inilarly oonatitutad area. I till another 
poaaibility tmuld ba to aatinata the 
p-vaotor diatribution for tha cxmtaxt olaa- 
■ifioation froM a "oonvantional* olaaaifiea- 
tien with "raaaonably good* aoeuraoy. Ml 
of thasa nathoda produea an aatinata of tha 
p-vaotor diatribution. and a'"ofttoial guao- 
tion on «dtioh hingaa tha utility of thin 
approach ia how nenaitiva the cmtaxtual 
algorithm ia likely to ba to tha *goodnaaa* 
of tha aatinata. H)ia quantion ia tha 
Bubjact of ongoing raaaaroh. 

An axparimant was formulated to obtain 
aorna avidenoa eonoerning tha faaaibility of 
applying tha oontaxt olaaoifiar to a raal 
data nituaticMi. Tha data eat uaad oovared 
a aonatdiat larger araa of Grand Rapida, 
Michigan, mntaining both data aata 2a and 
2b. Data from amall araaa of known ground 
oovar wara uaad to aatinata tha training 
olaaa atatiatioa. and data from a dia joint 
aat of araaa of known ground oovar were 
uaad an "tant aamplaa* to evaluata tha 
olaaaifiar accuracy (unfortunately, tha aat 
uaad for thia tant waa rathar amall, 
conaiating of only 138 pixala diatributad 
among 4 urban elaaaaa). 

A non-oontaxtual elaaalficatlon wan 
parformad and found, baaad on tha teat aat, 
to ba 81.6 percent accurate. Tha p-vaotor 
diatributiona wara eatimatad from thia 
olaaaifieation and uaad to perform contex- 
tual claaaificationa uning four and eight 
nearaat neighbors. The four-neighbor 
olaaaifieation wa.q 83.1 percent accurate; 
the eight-neighbor olaaaifieation waa 84.6 
percent accurate. For this case, then, 
some improvement in classification accuracy 
waa achieved by Incorporating context in 
tha decision process, although tha improve- 
ment was not as dramatic as for the simu- 
lated data sets. Whether thia ia due to 
poor estimation of the p-vector distribu- 
tions or simply to less contextual informa- 
tion in tha overall data set will be 
established by further investigation. 


III. MULTIPROCESSOR l.MPLEMENTATION 
OF CLASSIFICATION ALGORITHMS 

Classification algorithms such as the 
context classifier (and even much simpler 
algorithms used for remote sensing data 
analysis) typically require large amounts 
of memory and computation. These ar« said 
to be processor bound. Sines many avail- 
able systems, such as the IBM 360/370 VM or 
the PDP 11/70 UNIX, ara used on a time- 
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■haring baala# a large prooaaaor-bound 
program foroaa tha operating ayatem to 
operate with leaa overall memory, Coroing 
the memory management to awap larga amounta 
of information in and out of main memory, 
■n^ia reduoea the effieieney of piooeaaing, 
forcing the prooeaaor to take longer on all 
joba involved. For example, when UHZX ia 
under heavy loada (typically three to five 
proceaaor-bound joba with 35 to 40 on-line 
uatra), the cru apenda \xp to COI of ita 
time on operating ayaten taaka auch aa 
nwmory management. One way to apead up the 
proceaaea would to add a dedicated 
■pecial-purpoae proceaaor. Through tha uae 
of parallel prooeaaora, the ayatem through- 
put could be increaaed even more. Varioua 
dedicated ayatema have been propoaed, auch 
aa pipelined prooeaaora* , multimicrocomouter 
ayatema^'*, and apeoial purpoae ayatema^. 

To demonatrate the uae of a auch a 
ayatem on a teak leaa complex than tha 
contextual claaaifier, conaider the analyaia 
of Landaat data uaii 7 a Bayea maximum likeli- 
hood claaaifier (MU ) . tandaat meaaurcmenta 
are taken from four apectral banda and 
received aa a data vector. Baaed on doci- 
aion theory akin to that developed in the 
previouB aeotion the vector ia elaaaified by 
determining the probability that it belonga 
to each information claaa and aaaigning it 
to the claaa for which thia probability ia 
maximum. In thia caae, one approach would 
be to have one proceaaor compute the prob- 
ability for each of the claaaea. Such a 
method of proceaaing would yield a aubatan- 
tial increaae in throughput over a dedicated 
single-proceaaor ayatem. 

Conaider, for example, the Control Data 
Corporation (CDC) multiproceaaing ayatem 
conaiating of an array of dynamically micro- 
programmable proceaaora called Flexible 
Processors (FPs) *»**»**. The CDC PP cur- 
rently has no hardware facilities for 
floating-point operations, a disadvantage of 
the system. But the paralleliam of the 
system more than outweighs this fault. The 
basic clock cycle time is 125 nsec, but 
since the FP is designed to be connected to 
as many as 15 other FPa in a parallel and/or 
pipelined fashion, the effective throughput 
can be drastically increased, resulting in a 
potential effective cycle time of less than 
10 nsec. The CDC FP has been considered for 
its use in a large-scale image processing 
system’*. Its use in implementing a Bayes 
maximum likelihood classifier is demon- 
strated below. The techniques described are 
to be extended to the contextual classifica- 
tion algorithm. 

A configurational diagram of the FP is 
shown in Figure 3. (This is a preliminary 
FP design, but the final version should be 
very similar.) One of the features of the 


PP ia tha deiAla-bua arohitaetura whloh 
allows tha uaar to nanipulata data In IC-bit 
units, usa of IC-bit Intagar fonuta 
doublaa tha affaotiva storaga oapaoity or 
tha maohina, but 32-bit langths also ara 
aaaily handlad, tdiich makas it possibla to 
work with tha IBM 360/370 floating-point 
nunbara as wall aa tha POP 11/70 formats. 
Furthar, it is possibla to iag>laiMnt 
floating-point oparations in softwara, so 
tha Mohina is oapabla of doing floating- 
point arithnatie as is raguirad by tha elas- 
sification algorithms. 

in aach FP thara ara two ragistar files, 
one called tha temporary ragistar file and 
tha other tha larga register file. Both ara 
divided into 16-bit subunits. If tha needed 
path width is 16 bits, tha two files oan 
act like four files, thus creating more 
addraaaable uaar apace. A special feature 
of the temporary file is its separate read 
and write address registers, which can save 
much CPU time in many types of matrix 
oparations. it is possibla to do either a 
read or a write to either file and simulta- 
neously increment (or daeramant) tha address, 
further increasing throughput. Tha tempo- 
rary file is 16 words by 32 bits wide, while 
tha larga file is 4096 words by 32 bits wide. 
All of tha register files consist of 60 nsec 
random aoeasa memory (BAM) . 

Thara ara three general purpose regis- 
ters ( 6 PHs) called the E, F and G registers. 
All of these regiatera are connected to the 
arithmetic logic unit (ALU) • The E and G 
registers are readable only through the AbU. 
It is possible to shift the GPRa separately 
aa well aa combining the H and F regiatera 
to do a double-length shift.. The output of 
the ALU ia treated aa a register which is 
accessible in eight-bit units. Separate 
from the ALU ia a hardware integer multi- 
plier, which takes two eight-bit numbers and 
multiplies them to produce a 16-bit product. 
The input registers are the P and Q regis- 
ters, which are aach 16 bits wide. The 
user can choose which of the two groups of 
eight bits are to be multiplied. 

The FP is equipped with four index 
registers and four corresponding index- 
compare registers. There are four general 
compare registers called maintenance compare 
registers. All of the above are used for 
looping and can be incremented or decremen- 
ted during any statement not accessing 
those registers. 

The FP is equipped with a jump stack, 
so it is capable of handling standard types 
of program calls such as subroutine jumps. 
This stack is only 10 bits wide, and is, 
consequently, not suitable for storing 
data . 
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Xnput/oufeput (I/O) for th« fp dapcnda 
cm th« ovoraLl tyofcMi th« F7 array 

and ita hoafc naohina). Dlracfc I/O amonQ 
PPa and/or tha hoafc la dona via fcha M0» 

All I ARO and ARl ragiatara. niara ara 
intarlinkad naaiory unita on tha fA ayatan 
%d)ieh ara aooaaaad via tha X^n ’‘out 
ragiatara. intarrupta ara handiad tltrough 
tha Xntr (Xntarrupt) ragiatara, ao I/O ia 
fairly aaty and vary faat (about 3.1 awga- 
l^ftaa par aaoond) . 

Tha buaaaa ara connaetad to a ragiatar 
pair csallad tha BRO pair. Thaaa ara linkad 
to tha panal lighta on tha machina and can 
ba uaad for braakpointing or aa a OPR 
during axaoution. 

Pigura 4 ahowa note tha PPa ara linkad 
to tha boat and to aaoh othar. Tha ahift 
natteork ia ona waana of intar-PP communi- 
cation* tha cthar baing through intarlinkad 
manoriaa. Bach PP can addraaa certain 
mamory banka* tdtioh car. ba aooaaaad by 
certain othar PPa. Tha aharad iMimory (160 
naac cycle time) ia aapacially uaaful whan 
it ia naoaaaary to tranafar large aenounta 
of data between PPa. 

Pigura S ahowa a coding form for tha PP 
which ahowa* for axampla, that it ia poaaibla 
to conditionally inoramant an index ragiatar* 
do a program jump* multiply two aight-bit 
integara » and do a logical operation on the 
E and G regiatara* all aimultaneoualy. Thia 
type of operational overlap in conjunction 
with tha uaa of many prooeaaora executing 
concurrently greatly increaaea the effective 
apeed of the PP array. 

The ability to do a faat matrix multi- 
ple ia at tha heart of efficiently implemen- 
ting the Bayaa maximum likelihood claaaifier. 
The form for tha matrix multiplicationa ia‘ i 

(X-Uj^)’^ (C^^) (X-Uj^), 


covariance matrix to the data net aavaral 
would ba applied. Thia will* of oouraa* 
incraaaa tiia axaoution time by a factor of 
approximataly m/p. 


In atandard arithmatle* ona would firat 
multiply (X-U 4 )« and Cjl, creating a new 
vector. Thin vector trauld than be multi- 
plied by (X-Ui.) ranulHng in a aealar. in 
our impiMiantatlai* tha wdor haa bean 
aomawhat altarad. (X- 04 )C ia mulbipliad by 
a column of CT** aeoumuXating tha raaulta 
in a variable called "aum." After thia in 
dona for column j of Cjl* "aum* ia multi- 
plied by (X-Ui) 4 (the 3th element of (X-U|,»* 
accumulating thi raault in a variable 
called "hold" and ra-initialiaing "aum" to 
Tha following in a "pidgaon AtiOOb" 
deacription of tha prooeaa for ona pixel » 
total*0 

for j^l to n do 
begin I 
aum«0i 

for k-l to n do , 

aum>aum4DC k3 *C< [k * j ] j 
hold«holdtaum*DC j 3 1 

and I 


n ■ dimenaion of covariance 
matrix 

o[k3 " kth element of (X-Ui.) * 

. ccxnputed whan X ia loaded 

C Ck#j3 “ element in the kth row and 
jth column of 

At the end of the routine* tha value con- 
tained in tha "hold" variable ia the deairad 
aealar. Thia algorithm raguirea fewer 
atorea and fetchen than the atandard algo- 
rithm* no it ahortena the run time of the 
proceaa. All pointera are kept in the 
index regiater* further aimplifying the 
proceaa. Pinally* becauae only two 
accumulatora aro uaed* the three GPRa can 
be kept free for the floating-point opera- 
tiona* while the accumulatora are atored 
elaewhere. 


where X ia the data vector* ia the mean 
vector for the ith claaa* and C 4 ia the 
covariance matrix for the ith claaa. 

Conaider the uae of the PP array to 
perform theae olaaaificatlona. Aaaume 
there aro m diatinct claaaea and the computer 
ayatem containa p PPa. Bach PP ia aaaigned 
to proceaa m/p claaaea. The large file in 
each PP ia Initialiaed with the inverae of 
the covariance matricea and mean vectora for 
each claaa it waa aaaigned. The current 
data vector ia atored in each PP in the 
temporary file. When a new data vector ia 
loaded into an PP it overwritea the previoua 
one. Por aimplicity* but without loaa of 
generality, in the following aaaume that 
m " p. If m ia greater than p* then in each 
PP inatead of applying juat one inverae 


One way to perform thia algorithm ia 
to have the boat initially aend Cjl and U 4 
to PP i. The boat then aenda the current 
data vector X to PP 0* then to PP I* PP 2* 
etc. Au aoon aa the PP receivea the data 
vector* it begina the calculation of the 
value of the diacrimlnant function. After 
the boat given all PPa the data for pixel 
(i* j) * it waita until PP 0 haa calculated 
the value for ita diacrx.ninant function. 

The boat then retriovea the value of the 
diacrlminant function and loada PP 0 with 
the data vector for the next pixel. The 
host executea thia proceaa for all the PPa. 
When the laat PP haa tranamitted tha reault, 
the boat doea a compare and atorea the claaa 
index correaponding to the maximum of the 
diacrlminant valuea computed for thia pixel. 
Thua* the comparea are done by the boat 
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whll* fchft FP« ar« coMputlnq fch* diaoriminaat 
Tunofciona for th« naici: pixol, Minimiiino 


•tOt, it may ba poatibla to davalop a 
highly aeourata claaaifiar for oontaxt- 
rich aoanast 


Allowing 40 PP maohina oyelaa for eaoh 
floating point additicm and 9 FP maohina 
cyclaa for aaeh floating point multiply”, 
tha numbar of maohina oyolaa ia aa followai 

whara j » numbar of pixala and 
n m nvMbar of maaiurtMnti 
(aita of data vaetor) i 

aatup and elaar ragiatarat 9 

load maant 2n 

load oovarianea matrixt 4n’ 

load and normaliza data vactori 42jn+) 
inner loop of algorithm! 56jn^ 

outer loop of algorithm! 61jn 

56jn’ + 103jn + 4n* + 2n + j + 9. 

7hia aaaumaa that m, tha numbar of olaaaaa, 
aquaia p, tha number of prooaaaora. if m 
ia greater than p, tha runtima may ba 
approximated by multiplying by Im/pl, 

Exact compariaona of tha FP array with 
other ayatama are difficult without detailad 
Information about faotora auch aa pra*or 
poat-procaaaing done by tha boat machine 
and tha data praciaion uaed. However, to 
give a general idea of tha affactivaneaa 
of thia approach, conaider a 256 x 256 
claaaification of bandaat data (n«*4) uaing 
16 claaaea and a complate array of 16 FPa. 
Tha total procaaaing tima ia approximately 
10.7 aec. ESI.” atataa that their array 
proceaaor givea up to an incraaaa of 25 
timea over tha IBM 370/158. On tha claaai- 
fication of four channela into eight 
claaaea, their tima ia 6.3 aec. 

Due to the organization of the FP 
and tha fact the uaer can microprogram it, 
accurate mathematical analyaea of FP 
algorithms are complex. In order to study 
those timing questions, a simulator for a 
single FP was developed”. It has now been 
expanded to handle multiple FPa. The 
maximum likelihood olassifi'ir la currently 
being implemented on tho simulator to 
confirm the analytical timing results and 
to provide a working classifier program 
written in FP assembly language which 
could be run on the actual FP hardware. 

This will allow an accurate co8t-effective- 
neas study. 


IV. CONCLUSIONS 

The preliminary results from the use 
of context in classification are promising. 
By studying ways of estimating the p- vector, 
choosing the size and shape of neighborhood, 


The disoussion of performing of the 
neximiMi likelihood classifier demonstrates 
one way in which a multiple processor 
system can be used to speed up the 
processing of image data. The inq>lementa- 
tion of the classifier on the simulator 
and eventually on the eetuel F» eyst«n 
will provide hard data to verify the 
effectiveness of this approach. 

Through tha use of parallel, pipelined, 
and/or special purpose computer systems, 
such as the CDC Flexible Processor system, 
the types of computations required for the 
context classifier and other ocxnputation- 
ally demanding processes can be implemented 
efficiently. This will not only reduce 
the oong>utation time required to do contex- 
tual classification but will as well allow 
the investigation of techniques which may 
otherwise be considered infeasible. 
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Figure 1. A p-pix#l neighborhood. 
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Pigur® 3. Data path organisation in th® CDC Pl®xible Proc®®ior 
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Figura 4. FF array block diagram". 



Fiqure 5. Flexible Proceaaor coding form". 
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